SPhT-T07/028 



The centrality dependence of elliptic flow, the hydrodynamic limit, and the viscosity 

of hot QCD 

Hans- Joachim Drescher, 1 Adrian Dumitru, 2 Clement Gombeaud, 3 and Jean- Yves Ollitrault 3 

1 Frankfurt Institute for Advanced Studies (FIAS), Johann Wolfgang Goethe- Universitdt, 
Max-von-Laue-Str. 1, 60438 Frankfurt am Main, Germany 
2 Institut fur Theoretische Physik, Johann Wolfgang Goethe- Universitdt, 
Max-von-Laue-Str. 1, 60438 Frankfurt am Main, Germany 
'''Service de Physique Theorique, CEA/DSM/SPhT, CNRS/MPPU/URA2306 
CEA Saclay, F-91191 Gif-sur-Yvette Cedex. 
(Dated: February 1, 2008) 

We show that the centrality and system-size dependence of elliptic flow measured at RHIC are 
fully described by a simple model based on eccentricity scaling and incomplete thermalization. We 
argue that the elliptic flow is at least 25% below the (ideal) "hydrodynamic limit", even for the 
most central Au-Au collisions. This lack of perfect equilibration allows for estimates of the effective 
parton cross section in the Quark- Gluon Plasma and of its viscosity to entropy density ratio. We 
also show how the initial conditions affect the transport coefficients and thermodynamic quantities 
extracted from the data, in particular the viscosity and the speed of sound. 
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When two ultrarelativistic nuclei collide at non-zero 
impact parameter, their overlap area in the transverse 
plane has a short axis, parallel to the impact parameter, 
and a long axis perpendicular to it. This almond shape 
of the initial profile is converted by the pressure gradient 
into a momentum asymmetry-so that more particles are 
emitted along the short axis [l| . The magnitude of this 
effect is characterized by elliptic flow, defined as 



v 2 = (cos 2^ - $ R )), 



(1) 



where if is the azimuthal angle of an outgoing particle, 
$7? is the azimuthal angle of the impact parameter, and 
angular brackets denote an average over many particles 
and many events. The unexpected large magnitude of 
elliptic flow at RHIC [2| has generated a lot of activity 
in recent years. 

Elliptic flow results from the interactions between the 
produced particles, and can be used to probe local ther- 
modynamic equilibrium. If the produced matter equili- 
brates, it behaves as an ideal fluid. Hydrodynamics pre- 
dicts that at a given energy, vi scales like the eccentricity 
e of the almond P, Q . It is independent of its transverse 
size R, as a consequence of the scale invariance of ideal- 
fluid dynamics. If, on the other hand, equilibration is 
incomplete, then eccentricity scaling is broken and V2/S 
also depends on the Knudsen number K = X/R, where 
A is the length scale over which a parton is deflected by 
a large angle. 

Here, we show that the centrality dependence of V2/S, 
for both Au+Au and Cu+Cu collisions, can be described 
by the following simple formula [J|: 

hydro -, 

e ' e 1 + K/Kq ' 1 ' 

vijs is largest in the hydrodynamic limit K — ► 0. The 
first order corrections to this limit, corresponding to vis- 
cous effects, are linear in K. For large mean-free path, 



far from the hydrodynamic limit, V2/S ~ 1/K vanishes 
like the number of collisions per particle. One expects 
the transition between these two regimes to occur when 
A ~ R, hence that K ~ 1. A recent transport cal- 
culation [j| in two spatial dimensions indeed obtained 
K ~ 0.7.' 

Elliptic flow develops gradually during the early stages 
of the collision. Due to the strong longitudinal expansion, 
the thermodynamic properties of the medium depend on 
the time r, of course. The average particle density, for 
instance, decreases like 1 /r (if their number is approxi- 
mately conserved, see recent discussion in [6|): 



P(r) 



1 dN 
tS dy ' 



(3) 



where dN/ dy denotes the total (charged + neutral) mul- 
tiplicity per unit rapidity, and S is the transverse overlap 
area between the two nuclei. The quantities that we shall 
extract from V2 should be intepreted as averages over the 
transverse area S, and over some time interval around 
R/c s , which is the typical time scale for the build-up of 
V2 in hydrodynamics 4] . c s denotes the velocity of sound. 

The Knudsen number K is defined by evaluating the 
mean free path A = l/ap (a is a partonic cross section) 
at r = R/c s . Thus, 



1 



a dN 



K S dy 



(4) 



The purpose of this Letter is to show that the centrality 
and system-size dependence of the data for V2 at RHIC 
is described very well by Eqs. © and This provides 
three important pieces of information. First, such a fit 
allows us to "measure" the Knudsen number correspond- 
ing to a given centrality, which quantifies how close the 
dense matter produced in heavy-ion collisions at RHIC 
is to perfect fluidity. Second, the extrapolation to K = 



2 



allows us to read off the limiting value for v 2 y m /s ex- 
tracted from the data; this is useful for constraining the 
equation of state (EoS) of QCD via hydrodynamic sim- 
ulations, and we shall also see that it exhibits a rather 
surprising dependence on the initial conditions. Finally, 
using Eq. (|4|), we can convert the Knudsen number into 
the typical partonic cross section a (and viscosity) in the 
Quark-Gluon Plasma (QGP). Since only the combina- 
tion KqgCs actually appears in Eq. ([2]), uncertainties in 
Kq or c s then translate into corresponding uncertainties 
of er. Unless mentioned otherwise, our standard choice 
is c s = 1/V3 ~ 0.58 (ideal quark-gluon plasma) and 
Kq = 0.7. Letting Kq = 1 and c 2 — 2/3[38j instead re- 
duces the estimated a by a factor of two; on the other 
hand, taking Kq = 0.5 and c 2 = 1/6 increases a by the 
same factor. 

For the elliptic flow, v%, we use PHOBOS data for 
Au-Au and Cu-Cu 8| collisions. The same analy- 
sis could be carried out using data from PHENIX [9( or 
STAR 10]. The initial eccentricity e and the transverse 
density (1/ S)(dN/dy) are evaluated using a model of the 
collision. Two such models will be compared. The re- 
maining parameters i^ 3 "^ and a are fit to the data. The 
first step is to plot V2/e versus (l/S)(dN/dy) [ll|. Such 
plots have already been obtained at SPS and RHIC 
and they are puzzling: while V2/S increases with cen- 
trality, it shows no hint of the saturation predicted by 
Eq. |(2]) for K/K ^ 1, suggesting that the system is far 
from equilibrium Q ■ On the other hand, the value of V2 
for central Au-Au collisions at RHIC is about as high as 
predicted by hydrodynamics, which is widely considered 
as key evidence that a "perfect liquid" has been created 
at RHIC [3. 

It was understood only recently that the eccentricity of 
the overlap zone has so far been underestimated, as the 
result of two effects. The first effect is fluctuations in ini- 
tial conditions 14]: the time scale of the nucleus- nucleus 
collision at RHIC is so short that each nucleus remains 
in a frozen configuration, with its nucleons distributed 
according to the nuclear wave function. Fluctuations in 
the nucleon positions result in fluctuations of the over- 
lap area. Their effect on elliptic flow was first pointed 
out in Ref. [ll. It was later realized by the PHOBOS 
collaboration^, [l6[ that the orientation of the almond 
may also fluctuate, so that in Eq. |T|) is no longer 
the direction of impact parameter, but the minor axis 
of the ellipse defined by the positions of the nucleons. 
These fluctuations explain both the large magnitude of 
V2 for small systems, such as Cu-Cu collisions, as well as 
the non-zero magnitude of V2 in central collisions, where 
the eccentricity would otherwise vanish. They have to 
be taken into account in order to observe the expected 
saturation of i^/e at high density mentioned above. 

The eccentricity is usually estimated from the distri- 
bution of participant nucleons in the transverse plane 
(Glauber model). More precisely, we assume here that 
the density distribution of produced particles is given 
by a fixed 80%: 20% superposition of participant and 



binary-collision scaling, respectively [l7| . For Au-Au col- 
lisions, this simple model reproduces the centrality de- 
pendence of the multiplicity reasonably well (we assume 
that charged particles are 2/3 of the total multiplicity, 
and that dN/drj ~ 0.8 dN/dy at midrapidity), while it 
underestimates it for central Cu-Cu collisions by about 
10%. 

At high energies a second effect which increases the ec- 
centricity is perturbative gluon saturation, which deter- 
mines the pj_-integrated multiplicity from weak-coupling 
QCD without additional models for soft particle produc- 
tion. High-density QCD (the "Color-Glass Condensate") 
predicts a different distribution of produced gluons, 
dN/d 2 r±dy, which gives a similar centrality dependence 
of the multiplicity (Tt} but a larger eccentricity [ID, [h| . 
When particle production is dominated by transverse mo- 
menta below the saturation scale of the denser nucleus, 
then dN/d 2 r±dy ~ min(rip al . t (rjj, np art (rjj) traces the 
participant density of the more dilute collision partner, 
rather than the average as in the Glauber model [lj|. 
Precise figures depend on how the saturation scale is 
defined [201 ]. Naively, the larger initial eccentricity pre- 
dicted by the gluon saturation approach is expected to 
require more dissipation in order to reproduce the same 
experimentally measured i>2. Somewhat surprisingly, we 
shall find that this expectation is incorrect, which under- 
scores the non-trivial role played by the initial conditions. 

Both effects, fluctuations and gluon saturation, were 
recently combined by Drescher and Nara [211 ]. In their 
approach, the saturation momenta and the unintegrated 
gluon distribution functions of the colliding nuclei are 
determined for each configuration individually. The fi- 
nite interaction range of the nucleons is also taken into 
account. Upon convolution of the projectile and target 
unintegrated gluon distribution functions and averaging 
over configurations, the model leads to a very good de- 
scription of the multiplicity for both Au-Au as well as 
Cu-Cu collisions over the entire available range of cen- 
tralizes. 

Having determined the density distributions of pro- 
duced particles from either model as described above, 
we obtain the eccentricity via [HI, H2] 

l~ 2 — ^K-^) 2 + 4< 

£ = y ( £ part) ) e part = 2~T 2 ' 

a x , o~ y are the respective root-mean-square widths of 
the density distributions, and a xy = xy — xy (a bar 
denotes a convolution with the density distribution for 
a given configuration while brackets stand for averages 
over configurations). The overlap area S is defined by 
S = 4na x a y We find it more appropriate to de- 
fine these moments via the number density distribution 
dN/d 2 r±dy rather than the energy density distribution 
dE±/d 2 r±dy. The reason is twofold: first, V2 is extracted 
experimentally from the azimuthal distribution of parti- 
cle number, not transverse energy; second, our CGC ap- 
proach describes the centrality dependence of the mea- 
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FIG. 1: Variation of the scaled elliptic flow with the density, 
assuming initial conditions from the Glauber model. The line 
is a 2-parameter fit using Eqs. ([2]) and Q. 
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FIG. 2: Same as Fig. [1] using CGC initial conditions. 



sured final-state multiplicity very well, which indicates 
that the ratio of final-state particles to initial-state gluons 
(including possible gluon multiplication processes [23]) is 
essentially constant. 

Figs. [1] and [2] display v 2 /e as a function of 
(l/S)(dN/dy) for Au-Au and Cu-Cu collisions at various 
centralities, within the Glauber and CGC approaches, 
respectively. For both types of initial conditions, Cu-Cu 
and Au-Au collisions at the same (1/ S)(dN/dy) give the 
same v^l e within error bars. Eccentricity fluctuations are 
crucial for this agreement [|[ . The figures also show that 
Eqs. @ and provide a good fit to the data, for both 
sets of initial conditions. On the other hand, the values 
of the fit parameters clearly depend on the initial condi- 
tions, which has important consequences for the physics. 

The first physical quantity extracted from the fit is the 



hydrodynamic limit, v 2 /s, obtained by extrapolating 
to (l/S)(dN/dy) -> oo. The values are v 2 ydm /e = 0.30 ± 



0.02 with the Glauber parameterization, and v 2 ydro /e — 
0.22±0.01 with CGC initial conditions. Comparing these 
numbers to the experimental data points one observes 
that deviations from ideal hydrodynamics are as large as 
30%, even for central Au-Au collisions. This is our first 
important result. 

So far, a quantitative extraction of the QCD EoS from 
RHIC data via hydrodynamic analysis was hampered 
by the fact that v 2 /e had not been factorized into the 
perfect-fluid part v 2 ydl ° /e and the dissipative correction 
1/(1 + K/Kq). For example, Huovinen found [26J that 
an EoS with a rapid cross-over rather than a strong first- 
order phase transition, as favored by lattice QCD [27l |. 
overpredicted the flow data. This finding was rather puz- 
zling, too, as it was widely believed that the RHIC data 
fully saturates the hydrodynamic limit. Our results sug- 
gest that ideal hydrodynamics should in fact overpredict 
the measured flow. That is, that one should not choose 
an EoS in perfect-fluid simulations that fits the data. 
Rather, the EoS could be extracted by comparing ideal 
hydrodynamics to v 2 ydm /e. 

The next result is that CGC initial conditions, which 
predict a higher initial eccentricity e, naturally lead to a 
lower hydrodynamic limit v 2 ydro /e. Now, close to the 
ideal-gas limit (c s = l/y/3), v 2 ydlo /e scales approxi- 
mately like the sound velocity c s Q. This means that 
CGC initial conditions imply a lower average speed of 
sound (softer equation of state) than Glauber initial con- 
ditions, by a factor of 0.22/0.3 ~ 0.73. 

The second fit parameter is the partonic cross section 
a. The larger a, the faster the saturation of v 2 /e as 
a function of (1/ S){dN/dy). For our standard values 
of Kq and c s we obtain a = 4.3 ± 0.6 mb for Glauber 
initial conditions and a — 5.5 ± 0.5 mb for CGC initial 
conditions. These values are significantly smaller than 
those found in previous transport calculations [241 ], but 
match the findings of ref. [25| . 

CGC initial conditions imply a larger value of a than 
Glauber initial conditions, that is, a lower viscosity. This 
can be easily understood. As already mentioned above, 
the CGC predicts a larger eccentricity e than the Glauber 
model for semi-central collisions of large nuclei (when 
there is a large asymmetry in the local saturation scales of 
the collision partners, along a path in impact-parameter 
direction away from the origin [l9l]). However, for very 
peripheral collisions or small nuclei, there is of course 
very little asymmetry in the saturation scales, and the 
eccentricity approaches the same value as in the Glauber 
model. This has been checked numerically in fig. 7 of 
ref. plj . and can also be clearly seen by comparing our 
figures: while in Fig.[5]w2/e for semi-central Au+Au col- 
lisions is lower than in Fig. [1] there is no visible differ- 
ence for peripheral Cu+Cu collisions. In all, with CGC 
initial conditions the scaled flow grows less rapidly with 
the transverse density, which is the reason for the larger 
elementary cross-section. 

The dependence of a on the initial conditions is prob- 
ably even stronger than the numerical values above sug- 
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gest, for the following reason. As alluded to above, our fit 
to the data really determines the product K o~c s , rather 
than a alone, ft appears reasonable to assume that Kq 
does not depend on the initial conditions. However, for 
consistency, the speed of sound c s entering the Knudsen 
number should match the one underlying the hydrody- 
namic limit v 2 ydro /e; hence, if CGC initial conditions re- 
quire a smaller c s by a factor 0.73, the elementary cross- 
section obtained above should be rescaled accordingly. 
This leads to our final estimate (Jcgc — 7.6 ± 0.7 mb. 

Our numerical results for a should be taken as rough 
estimates rather than precise figures, because of the un- 
certainties related to the precise values of K and c s . It 
is, however, tempting to convert them into estimates of 
the shear viscosity rj, which has been of great interest 
lately. A universal lower bound rj/s > 1/4tt (where s is 
the entropy density) has been conjectured using a corre- 
spondence with black- hole physics [28|, and it has been 
argued that the viscosity of QCD might be close to the 
lower bound. Extrapolations of perturbative estimates 
to temperatures T ~ 200 MeV, on the other hand, sug- 
gest that the viscosity of QCD could be much larger [2!| • 
On the microscopic side, rj is related to the scattering 
cross-section a. Following Teaney [3(J, the relation for 
a classical gas of massless particles with isotropic differ- 
ential cross sections (which applies, for example, to a 
Boltzmann-transport model) is r] = 1.264 T/a [31(. On 
the other hand, the entropy density of a classical ultra- 
rclativistic cas is s = 4n, with n the particle density, so 
that 



77 T 

- = 0.316 — 



0.316—. 

can c 



(6) 



The relevant particle density in Au-Au collisions at 
RHIC, which is estimated at the time when v 2 devel- 
ops [4(, is 3.9 fm -3 for both Glauber and CGC initial con- 
ditions, and T ~ 200 MeV. Our two estimates a — 4.3 mb 
(Glauber initial conditions) and a — 7.6 mb (CGC initial 
conditions) thus translate into A = 0.60 fm, rj/s = 0.19 
and A = 0.34 fm, rj/s = 0.11, respectively. These values 
for rj/s agree with those from ref. [32j if the mean- free 
path is scaled to our result, and also with estimates of 
rj/s based on the observed energy loss and elliptic flow 
of heavy quarks [33| , on transverse momentum correla- 
tions [34], or bounds on entropy production 0. Hence, 
for our best fit(s) rj/s is slightly larger than the conjec- 
tured lower bound, but significantly smaller than extrap- 
olations from perturbative estimates. On the other hand, 



our lower value is close to a recent lattice estimate [35 1 
for SU(3) gluodynamics, which gives rj/s = 0.134 ±0.033 
atT = 1.65T c . 

A complementary approach to incorporate corrections 
from the ideal-fluid limit is viscous relativistic hydrody- 
namics. A formulation that is suitable for applications to 
high-energy heavy-ion collisions has been developped in 
recent years [IfJ . A first calculation of elliptic flow [371 ] 



shows that for Glauber initial conditions and rj/s — 0.16, 
v 2 reaches about 70% of the ideal-fluid value for semi- 
central Au-Au collisions. It is interesting to note that our 
simple estimates are in good agreement with this finding. 
Using Eq. ([6]), rj/s = 0.16 corresponds to a = 5.1 mb, for 
which Eqs. © and (gj) give v 2 /v 2 ydro = 0.68. The com- 
parison to experimental data in Ref. [371 ] . however, ap- 
pears to favor lower values of rj/ s because the EoS used 
there underpredicts v 2 ydro /s ~ 0.3 required for Glauber 
initial conditions. Alternatively, simulations could be 
performed with CGC initial conditions which require 
only 4 ydro /e ~ 0.22. 

In summary, we have shown that the centrality and 
system-size dependence of the measured v 2 can be un- 
derstood as follows: v 2 scales like the initial eccentric- 
ity e (as predicted by hydrodynamics), multiplied by a 
correction factor due to off-equilibrium (i.e., viscous) ef- 
fects. This correction involves the multiplicity density 
in the overlap area, (l/S)(dN/dy). Two types of initial 
conditions have been compared: a Glauber-type model, 
and a Color-Glass Condensate approach. PHOBOS data 
can be described with both. In particular, there is good 
agreement between Cu-Cu and Au-Au data. The result- 
ing estimates for thermodynamic quantities and trans- 
port coefficients, on the other hand, depend significantly 
on the initial conditions. 

Color glass condensate-type initial conditions require 
lower viscosity and a softer equation of state (smaller 
speed of sound). The scaled flow extrapolated to van- 
ishing mean-free path is lower than for Glauber initial 
conditions by a factor of ~ 0.22/0.3 = 0.73; the effective 
speed of sound should also be lower by about the same 
factor. Our estimates for the viscosity are rj/s ~ 0.19 for 
Glauber initial conditions, and rj/s ~ 0.11 for CGC ini- 
tial conditions, but these numbers should be taken only 
as rough estimates. 

We have also shown that the data for the scaled flow 
indeed saturate at high densities to a hydrodynamic 
limit. In central Au-Au collisions at RHIC, v 2 reaches 
70% (resp. 75%) of the hydrodynamic limit for Glauber 
(CGC) initial conditions. The corrections to ideal hydro- 
dynamics are therefore significant, but reasonably small 
compared to unity, implying that (viscous) hydrodynam- 
ics should be a valid approach for understanding flow at 
RHIC. Also, the asymptotic limit of v 2 /e has been iso- 
lated and could now be used to test realistic equations of 
state from lattice-QCD with hydrodynamic simulations 
of heavy-ion collisions. 
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